VARIATIONAL SPLINES ON RIEMANNIAN MANIFOLDS WITH 
APPLICATIONS TO INTEGRAL GEOMETRY 



ISAAC PESENSON 



Abstract. We extend the classical theory of variational interpolating splines 
to the case of compact Riemannian manifolds. Our consideration includes 
in particular such problems as interpolation of a function by its values on a 
discrete set of points and interpolation by values of integrals over a family of 
submanifolds. 

The existence and uniqueness of interpolating variational spline on a Rie- 
mannian manifold is proven. Optimal properties of such splines are shown. 
The explicit formulas of variational splines in terms of the eigen functions of 
Laplace-Beltrami operator are found. It is also shown that in the case of in- 
terpolation on discrete sets of points variational splines converge to a function 
in C norms on manifolds. 

Applications of these results to the hemispherical and Radon transforms 
on the unit sphere are given. 



1. Introduction and Main Results 

In the present paper we develop variational interpolating splines in the context of 
a Riemannian compact manifold with an emphasis on the so called average splines. 

One of the basic examples of such manifolds is the unit sphere S 2 . The analysis 
on the two dimensional sphere S* 2 found many applications in computerized to- 
mography, statistics, signal analysis, seismology, weather prediction, and computer 
vision. During last years the interpolation problem on S 2 , the problem of evaluating 
the Fourier coefficients of functions on the unit sphere and closely related problems 
about quadrature formulas on S 2 attracted interest of many mathematicians. 

The theory of variational splines on R d can be found in [2] and [7] . The inter- 
polation theory on S 2 was initiated in [13] and [16]. An approximation theory on 
S 2 along with many applications and an extensive list of references can be found 
in the recent monograph [4]. An approach to interpolation on manifolds, which is 
different from our, was developed in [3]. 

Our paper was motivated by the following problem which is of interest for integral 
geometry. 

Let M, dimM — d, be a Riemannian manifold and M v , v = 1, 2, N, is a family 
of submanifolds dimM v = <i„, < d v < d. Given a set of numbers v±, v-i, Vn we 
would like to find a function for which 

(1.1) f fdx — v Vl v — 1, 2, N. 



1991 Mathematics Subject Classification. 42C05; Secondary 41A17, 41A65, 43A85, 46C99 . 
Key words and phrases. Riemannian manifold, Laplace-Beltrami operator, variational splines, 
hemispherical transform, spherical Radon transform. 

1 



2 



ISAAC PESENSON 



Moreover, we are interested in a "least curved" function that satisfies the previous 
conditions. In other words, we seek a function that satisfies (1.1) and minimizes 
the functional 

||(l + A)*/ 2 u|| 

for appropriate real t where A is the Laplace-Beltrami operator on M. Note that 
in the case when the submanifold M v is a point the integral (1.1) is understood as 
a value of a function at this point. 

Our result is that if s is a solution of such variational problem then the dis- 
tribution (1 + A)'s should satisfy the following distributional pseudo-differential 
equation on M for any ip £ Cfi°{M), 

(1.2) f ifj(l + Afsdx = J2 a » [ <P dx > 

where coefficients a v G C depend just on s. 

This equation allows one to obtain the Fourier coefficients of the function s with 
respect to cigen functions of the Laplace-Beltrami operator A. 

From the very definition our solution s is an " interpolant" in the sense that it 
has a prescribed set of integrals. Moreover, we show that the function s is not just 
an interpolant but also an optimal approximation to the set of all functions / in 
the Sobolev space H t (M) that satisfy (1.1) and 

(1.3) ||(l + A)*/ 2 /||<tf, 

for appropriate K > 0. Namely, we show that s is the center of the convex and 
bounded set of all functions that satisfy (1.1) and (1.3). 

In the Section 4 we develop an approximation theory by interpolating splines in 
the case when all submanifolds M v are points. We consider convergence of such 
interpolants in two different cases. In a first case a set of points M v — x v ,v = 
1, 2, N, gets denser and in the second one the order of smoothness of splines goes 
to infinity. For both types of convergence we give the rates of convergence of such 
interpolants in Sobolev and uniform norms on manifolds. 

Our approximation theorem in the case when the order of smoothness of in- 
terpolants goes to infinity but the set of points {x„} is fixed, leads to a natural 
generalization of the classical Sampling Theorem in the sense that it allows a com- 
plete recovery of band limited functions, i.e. finite linear combinations of eigen 
functions of the Laplace-Beltrami operator. 

Note that some generalizations of a classical Sampling Theorem in the case 
of periodic functions and periodic splines on the unit circle can be found in [5], 
[12] and [14]. In our case the circle is replaced by a general compact Ricmann- 
ian manifold and trigonometric functions by eigen functions of the corresponding 
Laplace-Beltrami operator. The paper [1] contains a Sampling Theorem in the case 
of S 2 which is different from our. A result about uniform convergence of spherical 
splines to smooth functions on spheres was obtained in [4] with a different rate of 
convergence and for a different definition of splines on S 2 which requires a so called 
admissible set of knots. 

In the last section, 5 we apply our results to the hemispherical and to the Radon 
transforms on the unit sphere S d in R d+1 . Namely, we treat the inversion of 
the hemispherical transform as a specific interpolation problem, where informa- 
tion about a function is a set of values of integrals over hemispheres from a finite 
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collection of hemispheres. A smooth odd function / on S d will be obtained as a 
uniform limit of a sequence of " interpolants" of / which have the same integrals 
over a family of hemispheres as the / does. 

The same approach we apply to the spherical Radon transform: we obtain a 
smooth function / on S d as the limit of a sequence of "interpolants" of / whose 
integrals over subspheres from a family of subspheres coincide with the integrals of 
/. In both cases our interpolants provide an optimal approximations to unknown 
functions. 

Let us give a more detailed account of main results. We consider a compact 
Riemannian manifold M and the corresponding Laplace-Beltrami operator A. It 
is known that the operator A is elliptic, positive definite and selfadjoint in the 
space L2(M) constructed using a Riemannian density dx. If M has a non-empty 
boundary we assume the Dirichlet boundary conditions. 

The Sobolev space H t (M),t e R can be introduced as the domain of the operator 
(1 + A)*/ 2 with the graph norm 



It is known that in this case the operator A has a discrete spectrum = Ao < 
Ai < A2 < and one can choose corresponding eigen functions tpo,<pi,... which 
form an orthonormal basis of L 2 (M). A distribution / belongs to H t (M),t E M, if 
and only if 



This L 2 — inner product can be also considered as a pairing between H_ t {M) 
and H t (M) and in this sense every element of i/_ t (M) can be identified with a 
continuous functional on H t (M). 

Let F v , v = 1, 2, N, be a set of distributions from a certain H_ to (M),to > 0. 
Our main assumption about the family of functionals {F v } is that the functionals 
F v can be "separated" in the following sense: 

Independence Assumption. There are functions ft v G Cq°(M) such that 



where b vll is the Kronecker delta. 

Note, that this assumption implies in particular that the functionals F v are 
linearly independent. Indeed, if we have that for certain coefficients 71,72, ■•■,7jv 

N 



\\f\\ t =\\{l + A) t /*f\\,feH t {M). 




< 00 



where Fourier coefficients Cj(f) of / are given by 




5>^ = o, 



then for any 1 < [i < N 



N 



= ^7^(^)=7„- 
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The families of distributions that satisfy our condition include 

a) Finite families of 8 functional and their derivatives. 

b) Sets of integrals over submanifolds from a finite family of submanifolds of any 
codimension. 

Variational Problem 

Given a sequence of complex numbers v — {v v },u = 1,2,..., TV, and a t > t we 
consider the following variational problem: 

Find a function u from the space H t (M) which has the following properties: 

1) F v (u) =v v ,v = 1, 2, N, v = {«„}, 

2) u minimizes functional u — > ||(1 + A)*/ 2 u||. 

We show that the solution to Variational problem does exist and is unique for any 
t > tg even without any assumption. But we need the Independence Assumption 
in order to determine the Fourier coefficients of the solution. The solution to the 
Variational Problem will be called a spline and will be denoted as St(v). The set of 
all solutions for a fixed set of distributions F — {F„} and a fixed t will be denoted 
as S(F,t). 

Given a function / £ H t (M) we will say that the unique spline s from S(F,t) 
interpolates / on F if 

F v {f)=F v {s). 

Such spline will be denoted as s*(/). 

From the point of view of the classical theory of variational splines it would be 
more natural to consider minimization of the functional 

u — ^ \\& t/2 u\\. 

However, in the case of a general compact manifolds it is cascr to work with the 
operator 1 + A since this operator is invertible. 

Our main result concerning variational splines is the following. 

Theorem 1.1. If every functional F u ,u — 1,2, ...,7V, belongs to -ff_ to (M), if the 
Independence and Reality Assumptions are satisfied and if t > to + d/2, then for 
any given sequence v — {v v }, v = 1, 2, ...N, the following statements are equivalent: 

1) s t (v) is the solution to the Variational Problem; 

2) s t (v) satisfies the following equation in the sense of distributions 

N 

(1.4) (1 + Afstiv) = a u {s t {v))K, t > t + d/2, 

v=l 

where a\{st{v)), aN{st{v)) form a solution of the N x N system 

N 

(1.5) Pu^a u {s t {v)) = Vll ,n = 1, N, 
and 



(1.6) 



OO 

J=0 
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3) the Fourier series of St(v) is the following 

oo 
1=0 

where 

N 

Cj (**(«)) =< st(v),<Pj >= (1 + Aj-)~* a v {s t {v))F\A<Pj). 

The statement that (1.4) is satisfied in the sense of distributions means, that for 
any V £ H t (M),t > t , 

r N 

/ (1 + A)* s t (v)ipdx = V a v (s t (v))F v (<p). 

It is important to note that the system (1.5) is always solvable according to our 
uniqueness and existence result for the Variational Problem. 

It is also necessary to note that the series (1.6) is absolutely convergent if t > 
to + d/2. Indeed, since functional F v are continuous in the Sobolev space H to (M) 
we obtain that for any normalized eigen function ipj which corresponds to the eigen 
value Xj the following inequality holds true 

\F v ( Vj )\ < C(M, F)\\(l + A) to ^ 2 (pj\\ < C(M, F)(l + A,) 4 "/ 2 . 

So 

\F v {^)F^ 3 )\<C{M,F){l + \ )^, 

and 

|(1 + Xj )~ t F\X}Pj)F u (<pj ) | < C(M, F)(l + Xrf**- 1 ). 
It is known that the series 

3 

which defines the (— function of the Laplacc-Bcltrami operator, converges if r > 
d/2. This implies absolute convergence of (1.6) in the case t > to + d/2. 

We show that for a given function / e H t (M) its interpolating spline s t (f) 
has the following important property which means that it is always an optimal 
approximation in the sense of Golomb and Weinberger [6]. Namely, if Q(F, /, t, K) 
is the convex bounded and closed set of all functions g from H t (M) such that 

1) F v {g)=F v (f),v = 1,2,..., N, 
and 

2) || 5 || t <A-,forarealA->||* t (/)||t, 

then s t {f ) is the center of Q(F, f, t, K). This means that for any g £ Q(F, f, t, K) 

(1.7) \\s t (f) - g\\ t < \diamQ{F, f, t, K). 

To formulate our Approximation Theorem in the case when the distributions F v 
are Dirac distributions 5 Xu at points x u we need the notion of a p-lattice. 

We will say that a finite set of points X p = v = 1, 2, N, is a p-lattice, if 

1) the balls B(x v ,p/2) are disjoint, 

2) the balls B{x u ,p) form a cover of M. 
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The set of functionals F — {F„} associated with a p- lattice X p is the set of Dirac 
distributions 5 V = 5 Xu ,v = 1,2, ...,7V on Q°(M). 

We have the following result about convergence of splines in uniform spaces 
C k {M). 

Theorem 1.2. There exist a constant p(M) > 0, so that for any t > d/2 + k,d = 
dimM, there exists a constant C{M, t) such that for any p-lattice X p with p < p{M), 
and any smooth function f the following holds true 

(1.8) \\s2- d+ t{f){x)-f{x)\\ CHM) < (C(M,t)p 2 ) 2md \\(l + A) 2md+t f\\ 
for all m = 0,1, .... 

Moreover, if f is cu-band limited, i.e. belongs to the span of eigen functions 
whose eigen values are not greater than to, then 

(1.9) \\s 2 m d+t {f){x) f(x)\\ ck{M) < (1 + uf (C(M,t)p 2 (l + ,)) rd ||/||. 

The inequality (1.8) shows that convergence in C k (M) takes place when p goes 
to zero and the index 2 m d + t is fixed. 

The second inequality (1.9) shows that right-hand side goes to zero for a fixed 
p- lattice X p as long as 

p< (C(M)(l+oj))- 1/2 

and m goes to infinity. 

The second statement is, in fact, a generalization of the classical Sampling The- 
orem to the case of a general compact manifold. 

In the last Section, 5 we use above results for an approximate inversion of the 
hemispherical and Radon transforms on spheres. In these cases the functionals F v 
are integrals over hemispheres and subspheres respectively. In these situations by 
applying our Approximation Theorem on the " dual" sphere we show that if the set 
of hemispheres (resp. subspheres) gets "denser" than interpolants converge to the 
original function in C k {S d ) norms. 

Note that some of these results in the case of the spherical Radon transform were 
obtained in [11]. 

2. Variational splines on manifolds 

In this section we prove the existence and uniqueness of the interpolating varia- 
tional spline on a compact Ricmannian manifold. We also describe Fourier coeffi- 
cients of such splines with respect to an orthonormal system of eigen functions of 
the Laplace-Beltrami operator. 

We assume that F = = 1,...,N, is a family of distributions from a 

H_ to (M),t >0. 

Given a sequence of complex numbers {v v }, v = 1, 2, N, and a t > to we show 
that the corresponding Variational Problem does have a unique solution. 

Theorem 2.1. If every functional F v , v = 1,2, ...,N belongs to H_ to (M),t > 0, if 
t > to, then the Variational Problem does have the unique solution for any sequence 
of values (vi, V2, ...Vn). 

Proof. Consider the set V t °(F) C H t (M),t > to, of all functions from H t (M) such 
that for every 1 < v < N, F v (f) = 0. 

Given a sequence of complex numbers (v\, V2, vjv) the linear manifold 
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V t (F,Vl,...VN),t > t 

of all functions / from H t (M) such that F v {f) = v v ,v = 1,...,N, is a shift of the 
closed subspace V°(F), i.e. 

V t (F, Vl ,...,v N ) = V t (F)+g, 

where g is any function from H t (M) such that F v (g) = v v , v = 1, 2, N. 

Consider the orthogonal projection go of g £ H t (M) onto the space V®(F) with 
respect to the inner product in H t {M): 

</l,/2 >H t (M) = < (l + A)*/ 2 /l,(l + A)*/ 2 /2 >L 2 (M) = 

/ (l + A) t / 2 / 1 (l + A)*/2/ 2dx . 

It is clear that s t (v) = g — g £ V t (F, v\, vn) is the unique solution of the 
Variational Problem. Indeed, to show that s t (v) minimizes the functional 

U ^ + Ay/ 2 u\\ 

on the set V t (F,vi, ...,Vn) we note that any function from V t (F, v\, ...,Vn) can be 
written in the form s t {v) + h, where h £ V®(F). For such a function we have 

\\(l + Af 2 (a t (v) + h)\\ 2 = 
11(1 + A)*/ 2 s t (t,)|| 2 + 2 < s t (v),h> Ht{M) +||(1 + A)^ 2 h\\ 2 . 
Since s t (v) = g — go is orthogonal to V° (F) we obtain 

||(1 + A)*/ 2 ( St ( v ) + ah)\\ 2 = ||(1 + A)*/ 2 Si(u)|| 2 + M 2 ||(l + A)*/ 2 /i|| 2 , h £ V t °(F), 
that shows that the function s t (v) is the minimizer. □ 

Corollary 2.1. A function u £ H t (M) is a solution of the Variational Problem if 
and only if it is orthogonal to the subspace V®(F) and F v (u) = v v ,v = 1, 2, .... 

The next Theorem gives the characteristic property of splines. 

Theorem 2.2. If in addition to conditions of the Theorem 2.1 distributions Fi,..., Fjy 
satisfy Independence Assumption then a function s t (v) £ H t (M),t > t is a solu- 
tion of The Variational Problem if and only if it satisfies the following equation in 
the sense of distributions 

N 

(2.1) (l + A) t s t (v)=Y / ^t(v))'K. 

v=l 

In other words, for any smooth tp 

N 

< (1 + A)*s t (u), V >l 2 (m)= a v (s t (v))F v y>). 

v=l 
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Proof. We already know that every solution of The Variational Problem is orthog- 
onal to V t °(F) in the Hilbert space H t (M) i.e. for any h G V t °(F) 

(2.2) 0=<8t(v),h> Ht{M )= [ (l + A) 4 / 2 St (z;)(l + A)*/2/i. 

According to our Independence Assumption there exist a set of functions 
d = {& v },v = 1,2, ...,N, from Cg°(M) such that F^(jd u ) = 5^, where <5„ M is the 
Kronecker delta. Then for any ip G C§° (M) the function 

N 



belongs to V t °(F) and because of (2.2) 



JV JV 

Q=<s t {v),i>-Y^F v {i>)$ v>Ht{M) = / (l+A) t / 2 St (t;)(l + A)*/2( ?A _ Y,F V {$)4 V ) = 
Jjl + AfsAv) • 



In other words, 

JV 



/ (1 + A) t s t (v)ip = / (1 + A)*s t (i;)0„da;. 



If we set 



a v (8t(v),d)= [ (1 + Ays t (v)0 v dx, 
Jm 

we obtain that (1 + A)*s t (v) is a distribution of the form 

TV 

where 

FlW) = FJ$). 

So every solution of the variational problem is a solution of (2.1). 

Note, that if ( = is another Cq°(M) — family for which i^(C^) = then 
we have the identity 

TV 

which implies that 

a v (s t - a^(s t (w),C) = 
i.e. coefficients aj,(s t (w), $) = a„(s t (v)) are independent of the choice of the family 
of functions "3. 

Conversely, if u is a solution of (2.1) then since F v belongs to the space Ht (M), 
and t > to > 0, the Regularity Theorem for elliptic operator (1 + A)* implies that 
u G H- to+ 2t(M) C Ht{M) and for any h G V t °(F) 
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N 

<u,h > Ht{M) =< (l+A)*/ 2 ?i, (l+Af/ 2 h >=< (1+A)*u, h >= M"W») = 0, 

that shows that u is a the solution of The Variational Problem. □ 

As a consequence of the Theorem we obtain the fact that the set of all solutions 
of The Variational Problem is linear. In particular, every spline s t (v) G S(F, t) has 
the following representation through its values F v (s t (v)) =v v ,u= 1, N, on X: 

N 

(2.3) s t (tO = ^>X, 

where F u (s t (v)) — v v , and l" G S(F, t), v = 1, 2, iV, is so called Lagrangian spline 
that defined by conditions F tl (l 1 ') — S vfl , /J, = 1, 2, AT. 

To obtain another representation of splines we will need the solutions £7* , of the 
following distributional equations 

(2.4) (1 + AfEt = F~ v . 

To find E l v we note that in the sense of distributions 



(2-5) F v = Y,F v (Vi)Vi 

3=0 

which shows that 

oo 

(2.6) Et = j2^ + ^r t F^p])^- 

3=0 

According to the Theorem 2.2 every spline s t (v) is a solution of 

N 

(l + A) t s t (v) = J2<*Ast(v))'K, 

and along with (2.4) it gives 

(1 + AfsAv) = a v (s t (v)) ]T K&jJvj, 

v 3 

that implies the following representation 

N 

(2.7) s t (v)=J2a v (s t (v))El. 

Note that so far we have used just the assumption that t > t - To get more 
information about s t (v) we will need a stronger assumption that t > t + d/2. 

The next Theorem shows how to find an explicit connection between a sequence 
of values (v\,...,Vn) and corresponding sequence (ai(st(v)), otN(st{v))) in the 
case when t > to + d/2. We also assume that Independence Assumption is satisfied. 
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Theorem 2.3. If = Ao < A 2 < is the sequence of eigen values of A and 

ipo, ipi, ... is the corresponding sequence of orthonormal eigen functions, then for any 
spline s t (v) G S(F,t), such that F = {F„}, F v G H_ to (M),t> t +d/2,F v (s t (v)) = 
v v ,u = 1,...,N, the vector a(s t (v)) — (ai(s t (v)), Q!jv(st(w))) is the solution of 
the following N x N system 

N 

(2.8) P vli a v (s t (v)) = w„, /1 - 1, N, 

v=l 

where 

00 

(2-9) /V - ^(1 + Aj )~ t F~^X}Pj)F u (tpj ) . 

Proof. From (2.6) 

c j (Ei) = (l + X j )- t F~{^), 

or 

00 

^ = 2(1 + ^)-*^)^. 
3=0 

Combining this formula with (2.7) we obtain 

JV 

N 

^2a v (s t (v))P VIJ ,, 

where 

00 

fa? = ]>> + ^r'FA^F^j). 

3=0 

□ 

Note that according to the existence and uniqueness result the system (2.8) is 
always solvable. 

It was explained in the introduction that the series (2.9) is absolutely convergent 
iit > to + d/2. 

The following Theorem gives the Fourier coefficients of splines. 

Theorem 2.4. If the Independence Assumption is satisfied and t > to + d/2, then 
the Fourier coefficients of the spline s t (v) 6 S(F,t) are given by the following 
formulas 

N 

cA s t( v )) =< st(v),<pj >= (1 + Xj)-* a v { s M)FAjpj),j = 0, 1, ... 

where the vector (ai(s t (v)), ajv(st(u))) is the solution of the corresponding sys- 
tem (2.8). 
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Proof. If 

OO 

s t(v) = ^Cjistiv))^ 

then 

OO 

3=0 

and at the same time by (2.1) 

N 

Combining last two formulas with (2.5) we obtain 



The Theorem 2.4 is proved. 



3. Another extremal property of variational splines 



□ 



The goal of the section is to show that variational splines provide an approxima- 
tion which is optimal. 

Recall that for a given family of distributions F — {F^.v = 1,2,..., TV, from a 
_ff_ t() (M), i > 0, and a set of complex numbers v\, v 2 , vjv the notation V t (F; v\, ...vn) 
means the linear manifold of all functions / from H t (M), t > t 0} such that F v (f) = 
v v ,v = 1,2, ...,N. 

Lemma 3.1. For any g e V t (F;vi, ...vn) 

1/2 



\g\\t > (^^v^istiv))^ 



where s t (v) is the unique spline from V t (F; V\, ...vn). 

Proof. Let us note that the distance from zero to the subspace Vt(F, v\, vm) in 
the metric of the space H t (M) is exactly the Sobolev norm of the spline s t (v) £ 
V t (F,vi, ...,vn). This norm can be expressed in terms of the sequences v\,...,Vn 
and ai(s t (v)) 7 aAr(s t (v)). Indeed, 



ll*t(tOlk(M) =< (l+A)*/ 2 St ( u ),(l+A)*/ 2 St ( u ) > 1 /2 = < (l+Ay St (v),s t (v) >V*= 

N / N \ 1 / 2 

<^2a u (s t (v))F u ,s t (v) > 1 ' 2 = [^a u {s t {v))v v 

v=l \u=l ) 

The Lemma is proved. □ 
It other words the intersection 

Q(F, f, t, K) = V t (F, vi, .... v N ) p| B t (K), 
where B t {K) is the ball in H t (M) of radius K, is not empty if and only if 
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K > \\s t (v)\\ Ht{M) = a v(st(v))v v 
\ v =i 

By the very definition spline s t (v) is a solutions of an optimization problem. 
Now we prove another extremal property of s t (v). 

Lemma 3.2. The function s t (v) is the center of the convex, closed and bounded 
set Q(F,v,t,K) for any K > \\s t (v)\\ Ht{M) . 

Proof. We will show that if 

s t (v) +h £ Q(F,v,t,K) 

for some function h from the Sobolev space H t (M) then the function s t (v) — h also 
belongs to the same intersection. Indeed the last assumption shows that for any v : 
F v {h) = and then by (2.2) 

f (l + A) t / 2 s t (w)(l + A) t /2^ = o. 

J M 

But then 

||(1 + A)^(s t (v) + h)\\ = ||(1 + A)^(s t (v)- h)\\. 
In other words, 

||(1 + A)^(s t (v) -h)\\< K 

and because F v {s t {v) + h) = F v (s t (v) — h) for any v = 1,2, ...,7V, the function 
s t (v) — h belongs to Q(F, v, t, K). 

The Lemma is proved. □ 

Corollary 3.1. For any g G Q(F,v,t, K) the following inequality holds true 
\\st{v) - g\\t < ^diamQ(F,v,t,K). 

4. Spline approximation on manifolds 

In this section we develop the approximation theory by variational splines in the 
case when the set of distributions Fi is a set of delta functions on certain set of 
points of M. 

Definition 1. We will say that a finite set of points X = (xi,...,Xn) is a p-lattice, 
if 

1) The balls B(xi,p/2) are disjoint. 

2) The balls B(xi,p) form a cover of M. 

We will need the following result from [8], [9]. 

Lemma 4.1. There exist constants C(M) > 0,p(M) > such that for any p < 
p(M), any p-lattice X p = {xi\ and for any f £ Fi2d{M) such that f{x v ) = for 
all x v £ X p , the following inequality holds true 

11/11 < C(M)p 2d \\{l + A) d f\\,d = dimM. 
The next goal is to extend the last estimate to higher Sobolev norms. 
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Theorem 4.2. There exist constants C(M) > 0,p(M) > 0, such that for any 
< p < p(M), any p-lattice X p = {x v }, any smooth f which is zero on X p and 
any t > 

\\(l + A) t f\\<(C(M)p 2d ) 2 ||(1 + A) 2 d+ 7ll,*>0 
for all m = 0,1, .... 

We will obtain this estimate as a consequence of the following Lemma which is 
similar to a Lemma from [10]. 

Lemma 4.3. If for some f G H 2s (M),a, s > 0, 

(4-1) H/ll <a||AV||, 

then for the same f, a, s and all t > 0, m = 2 , 1 = 0, 1, 
(4.2) HA'/II <a m ||A ms+t /||, 

iff G H^ ms+t ^{M). 

Proof. Let us remind the reader that {Xj} is the set of eigen values of the operator 
A and {<fj} is the set of corresponding orthonormal eigen functions. Let {cj =< 
f,<Pj >} be the set of Fourier coefficients of the function / with respect to the 
orthonormal basis {<fj}- Using the Plancherel Theorem we can write our assumption 
(4.1) in the form 

ll/ll 2 <« 2 ( E A fM 2 + E A fM 2 

Since for the first sum » 2 A 2s < 1, 

o< E (M 2 -« 2 AfM 2 )< E (« 2 AfM 2 -M 2 ). 

Multiplication of this inequality by a 2 A 2s will only improve the existing inequality 
and then using the Plancherel Theorem once again we will obtain 

11/11 <a||AV||<a 2 ||A 2 V||. 
It is now clear that using induction we can prove 

11/11 <a m \\A ms flm = 2 l ,leE. 
But then, using the same arguments we have for any t > 

< Yl ( fl2TA f "N 2 - a 2(m+T) A 2(m+T)s M 2 ) < 

A 3 <a-!/s 

(a 2(m+T) A 2(m+T)s |c,| 2 - a 2T X 2 / s \ Cj \ 2 ), 

A^a-Vs 

that gives the desired inequality (4.2) if t — sr. □ 

To prove the Theorem 4.2 it is enough to apply the last Lemma 4.3 to the Lemma 
4.1 with a = C(M)p 2d . 

Now we can formulate and prove our Approximation Theorem. 
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Theorem 4.4. There exist constants C{M),p{M) > such that for any < 
p < p{M), any p-lattice M p , any smooth function f and any t > the following 
inequality holds true 

||(1 + A)*( S2md+t (/) - /)|| < (C(M)p*) 2md ||(1 + Af m+t fl 

for any m = 0, 1, .... In particular, if f is an uj-band limited function, i.e. f is a 
linear combination of orthonormal eigen functions whose corresponding eigen values 
belong to the interval [0,w], then 

||(1 + A)\s 2md+t (f) - f)\\ < (1 + uf (C(M)p\l + u)f ll/H, 

where m = 0, 1, .... 

Moreover, if t > d/2 + k then there exists a C(M,t) such that 

\\(s 2 m d+t (f)(x) f(x))\\ cHM) < (C(M 1 t)p 2 ) 2md ||(1 + A) 2md+t /||,m = 0, 1, ... 
and respectively, 

\\(s2m d+t (f)(x) - f(x))\\ cHM) < (1 + uf (C(M, t)pHl + uj)) 2 11/11, m = 0, 1, 
if f is an u-band limited function. 

Proof. To prove first two inequalities it is enough to use Theorem 4.2, the mini- 
mization property of splines and the inequality 

||(1 + A)V||< (l + wWH 
for all / which are linear combinations of eigen functions whose eigen values are 
not greater uj. 

To prove the inequalities in the uniform norms we use the Sobolev Embedding 
Theorem. The Approximation Theorem 4.4 is proved. □ 

5. Applications 

Example 1 

In this section we illustrate our results in the case of the hemispherical transform 
[15]. 

We consider the unit sphere S d C R d+1 and the corresponding space L 2 (5" i ) 
constructed with respect to normalized and rotation-invariant measure. Every vec- 
tor £ <E S d defines a hemisphere C S d as the set of all vectors x € S d for which 
£ • x > 0, where £ • x is the standard inner product in R d+1 . The correspondence 

will be treated as a correspondence between the set of all hemispheres of S d and 
points of the dual sphere S d . 

We denote by YJ an orthonormal basis of spherical harmonics in the space 
L 2 (S d ), where j — 0, 1, i = 1,2, n d (j) and 

is the dimension of the subspace of spherical harmonics of degree j . 
The Fourier decomposition of / e L 2 (S d ) is 
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where 



iAf)= I fYjdx=<f,Yj> L2{sd) . 
J s d 



To every function / e L 2 (S d ) the hemispherical transform T assigns a function 
Tf € L 2 {S d ) on the dual sphere S d which is given by the formula 

= / f(x)dx. 

J£-x>0 

For every function / e L 2 (S d ) that has Fourier coefficients Cij(f) the hemi- 
spherical transform can be given explicitly by the formula 

T.f(0 = nW^mjCijWYjiCZ E S*, 
where mj =0, if j is even and 

r(i/2) 



(_!)(j-i)/2. 



r((j + d + i)/2))' 

if j is odd. 

The transformation T is one to one on the subspace of odd functions (i.e. 
f(x) = —f(—x)) of a Sobolev space H^ dd {S d ) and maps it continuously onto 

TTodd ( ad\ 

I1 t+(d+l)/2\ D * h 

T(Hf d ){S d ) = H° d + d d+l)/2 (St). 

Let {h u },i> = 1,2, ...,N, be a finite set of hemispheres on S d . We consider 
functionals F u on L2(S d ) which are given by formulas 

F K =F V = j fdx. 

We will assume that the set of points S = on the dual sphere S d that 

corresponds to the set of hemispheres h v is symmetric in the sense that S = — S. 
Under this assumption we choose a t > and an odd function / and consider the 
following variational problem: find a function s t (f ) € H t {S d ),t > such that 

1) F v (s t (f))=F v (f),u = 1,2,..., N, 

2) st(f) minimizes norm ||(1 + A)*/ 2 s t (/)||. 

Since S = — S and function / is odd, the solution s*(/) will be an odd function. 
According to the Theorem 1.1 the Fourier series of St(f) is 

st(f)=J2^(st(f))Yj, 

i,3 

where the Fourier coefficients Ci,j(st{f)) OI s t(f) are given by formulas 

N 

c i Mf))=<s t (f),Yl>=(l + X iJ )- t J2Mst(f)) / Yjdx, 
where vector a(s t (f)) is the solution of the following N x TV system 

N 

$>«W a t(/)) = / fdx,n = 1,2,. ..N, 
b vll = V(l + A^)-* / Yjdx f Yjdx. 

Jh v Jh u 



where 
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This spline provides the optimal approximation to / in the sense that it is the 
center of the convex set Q(F, f, t, K) of all functions tp from H t {S d ) that satisfy 

(5.1) [ i[>dx= [ fdx 



and the inequality 

(5.2) ||(1 + A)'/ 2 V>|| <K 

for any fixed K that satisfies the inequality 



N 

K>\\8 t {f)\\ t = Y i <*M{f)) / /<*» 



Our results about the hemispherical transform arc summarized in the following 
theorem. 

Theorem 5.1. For a given symmetric set H = {h u }, v = 1,2, N, of hemispheres 
h v , an odd function f and any t > define the function s t (f) by the formula 



where 



and 



s t (f)=J2*Ast(f))Yj, 

N 

CiAs t {f)) = {l + \ i j)- t Y l <*Mf)) / Y j 



3 dx, 



N 



'%2b vli a v (st{f))=v ll ,v li = / fdx,n = l,2,...,N, 

6^ = ^(1 + ^)-*/ Y]dx f Yjdx. 

The function s t (/) is odd and it has the following properties. 

1) Integrals of the function st(f) over hemispheres h v have prescribed values v v : 

(T« t (/)) (£„) = f s t {f)dx = v v ,v = 1,2,..., N. 

2) Among all functions that satisfy (5.1) function s t {f) minimizes the Sobolev 
norm 

/ N \ V2 

||(1 + A)*/V/)||= [Y,a v {s t {f))vA . 

3) Function St(/) is the center of the convex set Q{F, f, t, K) of all functions g from 
H t {S d ) that satisfy (5.1) and the inequality 

(5.3) IKl + A)*/ 2 ^ <K, 

1 /2 

for any fixed K > a u (s t (f))vjj . In other words for any g £ Q(F, f, t, K) 

\\s t (f)-g\\ t <^diamQ{F,f,t,K). 
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Our next goal is to estimate the rate of convergence of st(f) to / in situations 
when either the set S gets denser or the smoothness t of splines s t (f) gets larger. 

Note that we can also interpolate the hemispherical transform Tf on the set 
S = by constructing an odd spline 

St = S T (T(f)),T = t+(d+l)/2, 

which is the solution to the following minimization problem: 

1) Sr(&,) = T/(£„) 

2) s T minimizes the functional u — > ||(1 + A) T / 2 u||. 
The Fourier series of s T is given by the formula 

S T = y ' c i,j(Sr)Yj , 

where Fourier coefficients Cij(s T ) of s T are given by the formulas 

N 

Ci,j(s T ) =< s T ,Yj >= (l + \ij)- T Y,Oi v (sT)Yj{^), 

v=l 

where a{s T ) is the corresponding "jump" vector which is the solution of the follow- 
ing N x N system 

N 

J2 b^a„(s T ) = (T/)(^), n = 1, 2, ...N, 

where 

6^ = ^(l + A i;j )-F/(^)F/(^). 

i,j 

The spline s T is the center of the convex set Q(S, Tf, t,K),t = t + (d + l)/2, 
where S is the family of delta functionals {5 V } = {S^„},^ v G H, and 

k = K\\T\\ Ht ^ HT . 

Since Q(5, Tf, t, K) is the image of Q(F, f, t, K) under the linear hemispherical 
transform T, we obtain that 

Ts t (f) = s T = s T (Tf). 

Applying our Approximation Theorem we obtain the following result about con- 
vergence of interpolants in the case of hemispherical transform. In this Theorem 
we use the following parameter p 

p = supinf dist(£,„, / /i, 

as a measure of the density of the set 5. 

Theorem 5.2. There exists a constant C such that for any m = 0,1,..., any 
k < 2 m d — d/2 and for any odd smooth f we have 

\\s 2 M.f) - /life < (C P 2 ) 2md \\(l + &Y /2 Tf\\,m = 0, 1, ...,r = 2 m d +(d+ l)/2, 

and if f is an odd spherical harmonic polynomial of (odd) order < uj then 

II W/) - /IU < (C P 2 (1 + u)? md {l + u) {d+1)/2 \\Tf\\,m = 0, 1, .... 
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The first inequality shows that for any odd smooth function / the interpolants 
S2 m d(/) of a fixed order 2 m d,m = 0,1, ...,k < 2 m d — d/2, converge to / in the 
uniform norm C k (M) as long as the parameter p goes to zero, i.e. the set S on the 
dual sphere gets denser. 

The second inequality in the Theorem shows, that interpolants converge to an 
odd harmonic polynomial of order uj for a fixed set of hemispheres S if Cp 2 (1 +u) < 
1 and as m goes to infinity. This statement is an analog of the Sampling Theorem 
for the hemispherical transform. 

Example 2 

We discuss the spherical Radon transform. It associates to a function / on S d 
its integrals over great subspheres: 

Rf(e ± ns d )= f fdx, 

Je ± ns d 

where 9 n S is the great subsphere of S d whose plane has normal 9. 

If a function / G L 2 (S d ) has Fourier coefficients Cij(f) then its Radon Transform 
is given by the formula 

R(f) = n-^T((d + l)/2) ]T r jCiJ (f)Yj. 
where Yf are the spherical harmonic polynomials and 

rj = (-iy/ 2 T((j + i)/2))/r((j + d)/2) 

if j is even and rj = if j is odd. Because the coefficients rj have asymptotics 
(— l)- J / 2 (j/2)( 1 ~ d )/ 2 as j goes to infinity we have that R is a continuous operator 
from the Sobolev space of even functions H^ ven (S d ) onto the space H%+\2-i)/2(S d )- 
Its inverse is a continuous operator from the space H^^_ l y 2 {S d ) onto the space 
H* ven (S d ). 

Let v = 1, 2, N, be a finite set of equatorial subspheres on S d of codi- 

mension one and distributions F v are given by formulas 

FAD = I fdx. 

By solving corresponding variational problem we can find a spline s t (f) € H t (S d ) 
such that 

F v (a t (f))=F v (f),v = 1,2,..., N, 
and s t (f) minimizes norm ||(1 + A)'/ 2 s t (/) || . 

Because we are interested in even functions on S d it is natural to have even 
splines. So we will assume that the set of points H = on the dual sphere that 
corresponds to the set of subspheres w v is even in the sense that S = — S. It is 
clear that this assumption will force our splines to be even functions. 

Our results about the spherical Radon transform are summarized in the following 
theorem. 

Theorem 5.3. For a given symmetric p- lattice W = {w v } of equatorial subspheres 
{wis}, v = 1,2, N, an even smooth function f and any t > d/2 + k define s t (f) 
by the formula 

St(f) = J2cMf))Yj, 
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where 



and 



CiAst(f)) = (l + X i , j )- t f2Mst(f)) I Yjdx, 

N 

^b vll a v {st{f)) = v^v,,. = / fdx,n=l,2,...,N, 

u=l 

&„„ = YJX + \o)~ l I Y I dx ! Y I dx - 
The function s t (f) is even and has the following properties. 

1) Integrals of s t (f) over subspheres w v have prescribed values v v 

\ s t (f)dx = v v ,v = 1,2, ...,N; 

J w u 

2) among all functions that satisfy 1) function s t (f) minimizes the Sobolev norm 

1/2 



(i + A)*/ 2 St (/)|| = (|>„M/)h) 



3) function s t (f) is the center of the convex set Q(F, f,t, K) of all functions g 
from H t {S d ) that satisfy 1) and the inequality 

||(1 + A)*/ 2 5 || <K, 

1 /2 

for any fixed K > (j2^ =1 ct v {s t {f))v^ . In other words for any g G Q(F, f, t, K) 

\\s t (f)-g\\ t <^diamQ(F,f,t,K). 

An approximation result similar to the Theorem 5.2 can also be formulated. 

Using our Approximation Theorem we obtain the following result about conver- 
gence of interpolants in the case of Radon transform. In this Theorem we use the 
following parameter p 

p = supinf dist^,^),^,^ eS,i/ / /i, 
as a measure of the density of the set S. 

Theorem 5.4. There exists a constant C such that for any m = 0,1,..., any 
k < 2 rn d — d/2 and for any even smooth f we have 

\\s2Mf) - f\\k < {C P 2 ? md \\(l + A)^ 2 i?/||,m = 0, 1, ...,r = 2 m d + (d - l)/2, 

and if f is an even spherical harmonic polynomial of (even) order < uj then 

\\s2m d (f) - f\\ k < {Cp 2 {l+Lu)f md {l+uj)^/ 2 \\Rflm = 0, 1,.... 

The first inequality shows that for any even smooth function / interpolants 
s 2 m d(f) of a fixed order 2 m d,m — 0,1, ...,k < 2 m d — d/2, converge to / in the 
uniform norm C k (M) as long as the parameter p goes to zero, i.e. the set S on the 
dual sphere gets denser. 

The second inequality in the Theorem shows, that interpolants converge in the 
uniform norm C k (M) to an even harmonic polynomial of order u> for a fixed set of 
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subspheres S if Cp 2 {l + uj) < 1 and m goes to infinity. This statement is an analog 
of the sampling Theorem for the spherical Radon transform. 
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